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Abstract 

The existence of proper weak solutions of the Dirichlet-Cauchy 
problem constituted by the Navier-Stokes-Fourier system which char- 
acterizes the incompressible homogeneous Newtonian fluids under ther- 
mal effects is studied. We call proper weak solutions such weak solu- 
tions that verify some local energy inequalities in analogy with the 
suitable weak solutions for the Navier-Stokes equations. Finally, we 
deal with some regularity for the temperature. 
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1 Introduction 

We study the existence and the regularity of weak solutions for the Dirichlet- 
Cauchy problem constituted by the Navier-Stokes-Fourier system which is 
one of the prominent mathematical problems for full thermodynamical sys- 
tems describing flows of incompressible fluids (see for example [HI [9l [23] and 
the references therein). The presence of the temperature dependent viscosity 
into the momentum equation and the Joule effect into the energy equation 
are the main contributions on the nonlinear behavior of the coupled system 
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of partial differential equations. Other problem which arises from fluid ther- 
momechanics is the Boussinesq approximation describing the dynamics in 
the planet mantle (see [22] and the references therein). In fact, it is a sim- 
pler coupled system which does not include the Joule effect. For instance, in 
[T2l [5l [T7] even the viscosity is considered constant. 

Although recently the C(]0, T[; C°'"(fi))- regularity for any weak solution 
is found to the non-stationary Stokes system in 3D [21], the non-uniqueness 
and smoothness of weak solutions are related with the non-uniqueness of 
Leray-Hopf's solutions to the first initial boundary value problem for the 
three-dimensional Navier-Stokes equations. For any class of weak solutions 
of the non- stationary N-S equations in three-dimensional spaces, it is known 
that a weak solution cannot be used as a test function in the weak variational 
formulation. Indeed, in the N-S-F system, the viscosity is a temperature 
dependent function and the temperature is a solution of a parabolic equation 
with L^-data due to the existence of the Joule effect. By these reasons, in the 
work [6[ [25] (and some references therein) the N-S-F system is constituted 
by momentum and total energy equations in order to the dissipative term 
in the internal energy equation is obtained by approximative methods using 
the weakly sequential lower semicontinuity of the norm. On the other hand, 
in the N-S system the study of the class of suitable weak solutions is being of 
interest since [TJ [201 [SSI [27], if homogeneous Dirichlet boundary conditions 
are enforced (see for instance [T^ fTB]). 

In this work, we prove the existence of weak solutions to the problem 
under study and that among the weak solutions at least the existence of 
one proper weak solution is guaranteed, i.e. such weak solutions that verify 
some local energy inequalities. Moreover, we show a higher integrability 
of the gradient of the velocity, using the reverse Holder inequality, which 
is used in the proof of existence of weak solutions in order to deal with 
the L}- dissipative term in the energy equation. Some regularity for the 
temperature appears as a direct consequence. We refer to [11] for the proof 
of the higher summability of the gradient of the velocity in the stationary 
case. The regularity of strong two-dimensional solutions was solved by the 
author in the paper [TU] if some smallness on the data is considered. The 
study of the partial regularity for the velocity of the fluid will be an ongoing 
future work. 

Let n C M" be a bounded open domain sufficiently regular and T > 0. 
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1) 



Let us consider the boundary-value problem of the N-S-F system: 

^tu - div{fi{9)Du) + (u ■ V)u = f - Vp in Qr := ^^x]0, T[; 

div u = in Qt; 

dtO - diY{k{e)Ve) + u-Ve = n{e)\Du\'^ in Qt; (2) 

"lt=o = "0, d\t=o = ^0 in (3) 
u = 0, ^ = on(9fix]0,T[, (4) 

where p denotes the pressure, fi the viscosity, 9 the temperature, u the veloc- 
ity of the fluid and Du = |(Vu + Vu-^), f denotes the given external body 
forces, k denotes the conductivity and it is assumed constant, k{9) = k. No- 
tice that the assumption of a constant conductivity is not necessary for the 
proof of existence of weak solutions. It is a sufficient condition for the proof 
of existence of the local energy inequalities to the temperature solution. For 
simplicity, the constant density is assumed equal to one, and we do not con- 
sider the existence of the external heat source, since the heating dissipative 
term is the main mathematical difficulty. The product of two tensors is given 
hj D : T = DijTij and the norm by = D : D. 

The initial conditions are given in ([3]). For the sake of clarity we found 
convenient that the boundary conditions which are given in (jl]) are assumed 
homogeneous Dirichlet conditions. 

The outline of the paper is as follows: in next section we establish the 
appropriate functional framework and we present the main results. The 
Section [3] is devoted to the proof of the existence result (Theorem 12. ip . In 
Section m we prove some regularity result (Proposition 12.11) . 

2 Assumptions and main results 

Here we assume that Q C M" is a bounded open set such that its boundary 
dQ G C^. In the framework of Lebesgue and Sobolev spaces, for 1 < g < oo, 
we introduce [18] 



3l'%n) = {ue Wl'\Q) : V ■ u = in fi} 

with norm 

II ■ ||l,cj,!^ = II V ■ llq.f^) 

and we set 
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For any set A, we write (m, v)a '■= uv whenever u E L'^{A) and v G L'^' (A), 
where q' = q/{q — 1) is the conjugate exponent to q, or simply (-, ■) whenever 
there exists no confusing at all, and we use the symbol (■, ■) to denote a generic 
duality pairing, not distinguished between scalar and vector fields. We denote 
by bold the vector spaces of vector-valued or tensor-valued functions. 

The following assertions on data are assumed as well as the following 
assumptions on the physical parameters appearing in the equations are es- 
tablished. 

(Al) f : ^ is given such that f G L'^^^+">\Qt) with eo > 0; 
(A2) /i : M — > M is a continuous function such that 

< /i# < < fi*, Vs G M; (5) 

(A3) uo G L2(fi), 6*0 G L\n) such that 

div Uo = in Q; ess inf 9o{x) > 0. (6) 

Definition 2.1 We say that the triple {u,p,6) is a weak solution to the 
Navier- Stokes- Fourier (N-S-F) problem (QP-Q' in Qt if 

ueU:= L°°(0, T; V{n)) n ^^(0, T; Jo''(^^)), P G L("+2)/"(gr), 
0eS:= L°°(0, T; L\Q)) n ^"(0, T; W^'\Q)), q < 

n + 1 

dtu eX := L^{0, T; W-^'2(fi)) n L("+2)/"(0, T; W-^'("+2)/n(^))^ 



£ g(n + l) 2(n + 2) 
and satisfies the variational formulation 

{dtVi,w)+ I (/i(0)Du : Dv + (u ■ V)u • v) = / (f-v + pdivv) 
JQt ^ ' JQt 

Vv G L°°(0, T; W^'°°(fi)), u|,^o = Uq m (7) 
{dtB, I (kVd - ^u) • V(pdxdt = [ fi{9)\Du\^cf>dxdt 

Jqt ^ ' Jqt 

V0 G L°°(0, T; iyo^'°°(fi)), ^1,^0 = 9^ m n. (8) 
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The embedding X ^ L("+2)/"(0, T; W-i'("+2)/"(n)) occurs for n = 2,3. 
For every u E U ^ L^*^"^^)/"((5t), the convective term verifies (u • V)u e 
L(n+2)/(n+i)^g^-) and Consequently (u ■ V)u ■ u ^ L^{Qt). Moreover, the 
advection term u ■ V6' ^ L}{Qt) ii 9 E £ and q < {n + 2)/{n + 1). Remark 
that 6 e8^ L'i^'^+^^I''{QT) and e L^(Qt) for £ > 1, i.e. g > 1 if n = 2 
and q > 15/14 ii n = 3 {q > 2ra(ra + 2)/((n + 4)(n + 1))). 

Definition 2.2 VFe say that a weak solution {u,p,9) of the N-S-F problem 
is proper in the following sense. The local energy inequality holds 



Jn 2 jQt=Qx]0,i[ 

<2 ij,{e)(pDu: {{u- a) (g)V(p) + \u - af(p{dt(p + u ■ V(p) + 
JQt JQt 

+2 / pyju • Vipdxdr + f ■ (u — 8L)(p'^dxdT, (9) 
JQt JQt 

for all ip G C^{Qt), a.e. t g]0,T[ and for any a G M", and two more local 
energy inequalities hold 

(VC+^V') {x, t)dx + CkJ^ j^^^^l^^dxdr < 
<[ y^C + THd,^ + kAij + u-V^)+ [ fi{e)\Dn\ \ f"^ (10) 

JQt JQt Is + J ' 



1-e 



/" (1 + ^)^"^ (^tV' + A;A^ + u ■ V^) da;dr, (11) 



for any ( > and < ^ < 1, and for all ip G C^{Qt) such that ''p >Q, a.e 
tG]0,T[. 

Remark 2.1 Two local energy inequalities are stated for the temperature 
because from ^^) if we take C ^ 0"*" we obtain 9 G L°°{0,T; Ll^^{Q)) and 
from / fli]) after some computations we get V9 G L1^^{Qt). 
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Theorem 2.1 (n = 2,3) Under the assumptions (A1)-(A3), the N-S-F prob- 
lem defined by has proper weak solutions. Moreover, Vu G 'Lfl^l'^^\QT) 
forO <e< min{(4-n)/(3n),l/(n + 2),eo} and 9 > zn Qt- 

Some interior regularity is proved. 

Proposition 2.1 (Interior regularity) Let {u,p,6) be a weak solution in 
accordance to Theorem \2.1\ Then 6 G W1^^^^^{Qt) ■ In particular, 6 G 

Moc {Qt}- 

Henceforth C will denote different positive constants depending on the 
data, but not on the unknown functions u, p or 0. 



3 Proof of Theorem 12.1 

The proof of Theorem 12. II is split in each subsection. In sections l3. 11 13.21 and 
13.31 for reader's convenience we delineate the main arguments concerning the 
existence of approximate solutions (for details, see [B]). The sections 13.51 
and l3.6l are new and they are the main contributions for the desired existence 
result. 



3.1 The Faedo-Galerkin method 

For I/, £ > fixed, there exists {u^'^-'^,p^'^, ^^'^'^}Ar_jv/eN being of the form 
[30l Caratheodory Theorem] 



N 

i=i 

M 

where {(w-', is a basis of WQ'^(r2) x Wl'^{Vl) with P > n, J-'e is the 

continuous functional such that maps u G WQ'^(f2) into p G W'^''^{Q) which is 
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the solution of the homogeneous Neumann problem for the Laplace equation 
(see, for instance, [13]) 

eAp{t) = div u(t) in Q 
Vp{t) ■ n = on on 

/ p{t)dx = 0, 



e\\pm2,2,n < C(fi)||u(t)||i,2,n, Vu(t) e wI^\Q) 

4Pmi,r,n < Cin,r)\\n{t)l,n, Vu(t) G Wl'\n) n (n) , r > 1. 
The functions c^'^^ = {c^'^'\ ■ ■ ■ , %'^^) and d^'^'^ = {d^'^ , ■ ■ ■ , d^^^) solve 



the following system of ordinary differential equations, for every M, G N, 



at 

-(^,K'^0,V-w^-) = (f,w^'), j = l,---,iV; (12) 

^(^^'*^w;^') - (0^'*^^,(u^'^O, V«;^') + A;(V^^'^^ Vw;^') = 
at 

= (M^"^''Ol^u^''?,t^^), J = 1,---,M, (13) 



under the initial conditions u^, ' given by the projections of uq and the 
mollification 6^ of 6q (after extending by zero outside f2), respectively, onto 
linear hulls of the base's vectors. Finally, A^jy is the divergenceless part of 
the Helmholtz-mollification decomposition 



which satisfies, a.e. t g]0,T[, 



d 



Myivi) := (xu) *uj -Vh 



with uj denoting a mollifier with support in a ball of radii z/. 




if dist(x,9fi) < 2z/ 

1 elsewhere. 



and h is the Helmholtz-mollification decomposition, that is. 



Ah = div[(xu) * cij] in f2 



Vh ■ n = on dQ 



hdx 



0. 
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Moreover, the standard estimates hold, independently of M, [121 fH] 



N,M 



sup II u 

ie[o,T] 

+^l|Vp^'^-^|llQ,<||uo|||^ + C||f||^,Q^; 



sup 11^^'^-^ 
te[o,T] 



2 ^k\\yQNM\\2 



2,n 



<r llfl^l|2 

2,Qt — Iro Il2,f7 



+ C(iV); 
< C(iV); 



L2{0,T) 



(14) 
(15) 

(16) 
(17) 



Hence, the initial value problem f|T2|) - (fT3|) has a global- in-time solution and 
passes to the limit as M tends to infinity (A^ fixed) by standard compactness 
arguments, resulting 

4-(u^, w^') - (MJu'') ® u^, Vw^) + (/i(^^)Z}u^, Dw^) - 
at 

-(^,(u^),V-w^-) = (f,w^-), j = l,---,iV, a.e. tG]0,T[; (18) 
{dtO'', 0) - (e^M.iu^), V0) + A;(V^^, V0) = 
= (MO|Du^r,0), V0GL2(O,T;iJi(fi)). (19) 

Moreover, the minimum principle holds, i.e. > a.e. in Qt- 

In order to pass to the limit on A^, when tends to infinity, the estimates 

(IT3]) - (IT71) are no more valid. Thus, we recall the additional estimates (for 

details, see [6]) 
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|oo,Li(n) 



< fx*\\\Du^\'%,Q^ + T\\9o\\i,n + \Qt\/2; 



(20) 



We^Wlr.. < ^(/.^iiiDu^nii.Q. + ii^olli,^.) X ii^^ii^^-.-f^^^^^ 

< C{u+^- + l]\\u^\\l^^,^^^+Cml^^; (22) 



\OtU ||2,w-i'2(n) 



N\ 



|i,iy-i.9(Q) 



< ||V^^^||,,Q, + C(z/)||^^||,,Q,+/x#|||Du 



N\2\ 



(23) 



for every exponent 1 < g < 2 — n/(n + 1) (cf. [3l H]). Moreover we have the 
strong convergences 

(24) 
(25) 

Then, the initial value problem f|T8|) - f|T9|) passes to the limit as tends 
to infinity {e fixed), concluding the below quasi-compressible approximative 
problem. 



Vp^ -^Vpe in L^(Qt); 
/x(^^)|L'u^|2 ^ /i(^,)|Du,|2 in L\Qt). 
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3.2 The quasi-compressible approximative problem 

From Section [STT] for each e > 0, there exists (u£,pe, 6^) in {L°°{0, T; L^(r2))n 
L2(0, T; W^'^(fi)))xL2(0, T; iy2,2(^))x^ such that dtu, G L2(0, T; W-i'2(fi)) 



and dt9s e L\0,T] W-^''^{n)), and it satisfies 



(StUe, v) + / Vue : V (g) Muiue)dx + / fi{6s)Due : Dvfia; = 
if7 in 

= / (f - Vj9,) ■ vrfx, Vv G Wl'\n), a.e. t g]0, r[; (26) 
Jn 

{dtOe, [ M^iUe) ■ Wee<pdz + k [ V^, ■ V0rf2 = 

JQt Jqt 

= [ /i(^.)l^u,p0dz, V</)GL-(0,T;iyoi'^'(fi)); (27) 



s Vpe- W(t)dx + / div u,0rfx = 0, V0 G W''\VL)- (28) 

u,(-,0) = uo, ^,(.,0) = ^o. 

In order to pass to the hmit, when e tends to zero {v fixed), the estimate 
( |22|) is no more vahd. To estimate Pe independently of e we choose v = Vt^^ as 
a test function in (l26l) . where rj^ is the solution of the following homogeneous 
Neumann problem for the Laplace equation (see, for instance, [T3] ) 



^r]e{t) = p^{t) - j pe{t)dx in VL 
Vrie{t) ■ n = on 



r]s{t)dx = 0, 



which satisfies 



WVeirnhn < CmWPemln, a.e. t g]0,T[. (29) 
After some technical computations, it results (for details, see [6]) 

WpsWIq^ < c{u), 

and consequently 

||a,u,||2_w_,,.(^) <C(z/). (30) 
9 



Note that i^{e,)\Du,\^ i^{e^)\Du^\'^ holds in L\Qt) (compare to 
and it is not required the strong convergence of the pressure Vpe — >■ Vp^ in 
L^(Qt) (compare to fl2^ ) since fl28l) holds for (p = p^- Then, the initial value 
problem fl26l) - fl27j) passes to the limit as e tends to zero (z/ fixed), concluding 
the below regularized problem. 

3.3 The regularized problem 



For each z/ G N, there exists {uu.pu, 6^) in W x L^(Qt) x S such that dtu^ G 
L2(0,T; W-i'2(n)) and dt9^ G L^{0,T;W-^''^{n)), and it satisfies 



Now in order to pass to the limit in flST]) - fl52]) as u tends to infinity, neither 
(!23|) nor (!30|) are valid. Following the argument of [6], we decompose the 
pressure pi, into pi, := Pi,^i + Pu,2 such that the two particular pressures, p^^i 
and pu^2-, belong to bounded sets of L*^"+^)/"'((5t) and L'^{Qt), respectively, 
independently of z/. For each t g]0,T[, p^^i is the unique solution to the 
problem 



for all G W'^'^{Q) such that V0 ■ n = on dfl, and define p^^2 '■ = Pu — Pu,i- 
Then Pi, 2 solves at each time level 




(31) 



(32) 



u,(-,0) = uo, 9,{-,0) = 9o. 




(p,,i(t),A0) = (u,®A^,(u,)(t),DV0) 




for all G Vr^'^(f2) such that V0 G WQ'^(fi), and the following estimate 
holds 



\\Pu,l\\in+2)/n,QT < C] 



Pu,2\\'2,Qt ^ ^• 
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Thus we conclude the following uniform estimates 

||^!^u,^||i,w-i''!(Q) < C; \\dtdu\\i,w-'^''^{n) < C, (34) 

for 1/i = n/[q{n + 1)] + n/[2{n + 2)] < 1 since q < {n + 2)/{n + 1) and 
n < 4. Then under compactness arguments (u,p, ^) satisfies the limit 
variational formulation ([7]). However, with the above estimates we only ob- 
tain 



Ou ■ V(t)dxdt + k V9- Vcpdxdt > {fi{9)\Du\ 
Jqt 



for all (j) G C^((5t) such that > [6]. So we will need to prove an additional 
estimate. 

3.4 Higher integrability of Vu^, 

First let us remark that (l3Tl) - fl32l) may be rewritten as 

= /(f -v+p^div v)cix, a.e. te[0,T] Vv e Wo'^(l]), (35) 
Jo. 

u|j^o = uo in Vt; 

+ j (kV9u-0uM,iu,)y\/(l)dx = j fxi9,)\Du,\^(l)dx, (36) 
a.e. t G [0, T] V0 G W^o '°°(fi), ^|t=o = in n. 
Let us prove the required estimate. 

Theorem 3.1 Let (u^,p^,^^) e U x L^''+^'^/''{Qt) x S verify the system f31\}- 
(ED then Vu, G L2e^')(QT) for < e < min{(4 - n)/{3n), l/{n + 2),eo} 
and the following estimate holds: 

||VU;,||2(l+e),Q(2„,K) < C (||VU;,||2,Q(2o,2R) + ||u,,||2(n+2)/n,Q(2o,2R) 

1 /2 \ 
+ ||f||2{l+eo),Q{2o,2R) + \M(n+2)/n,Q{zo,2R) ) ' (^7) 

/or any cylinder Q{zo,2R) := i?(a;o, 2i?) x]to — (2i?)^,to[cc Qt- 
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Proof. We denote the points of the space-time cyhnder hj z = {x, t) and 
employ a shorthand notation dz = dxdt. We write -f^^^ v := ^ Jq(^^ ^ 
whenever v G L^{Qt), where Un denotes the measure of the n-dimensional 
unit ball in R". For every zq = (xo,to) ^ Qt and R> small enough such 
that Q{zo, 2R) CC Qt, in order to prove the higher integrability of Vu,^ (cf. 
[H [m EHl Gehring Lemma]), it is sufficient to show the following reverse 
estimate 

Jq{zo,R) Jq{zo,2R) 

Bi f , ,2 , B2 f , ,0 , 

-n- Jq{zo,2R) -rt Jq{zo,2R) 

+B3-f \i\^dz + R-f |p,|("+2)/"ci2(38) 

Jq{zo,2R) Jq{zo,2R) 

for some 5 G [0, 1[ and positive constants Bi,B2,B3, independent of \iy,Pu 
and considering that, for n < 4, 

^ L^("+2)/n(g^)^ 2(n + 2) ^ 2. 

|^^|3/2 ^ ^4(n+2)/(3n)^g^)^ 4(n + 2) ^ 2. 

f GL2(^+-)(Qt), eo>0; 

Thus, we take < e < min{2/n, (4 — n)/(3n),l/(n + 2),eo} = min{(4 — 
n)/(3n),l/(n + 2),eo}. 

Adapting the argument used in j^, let ip G C^(Q(2;o, 2i?)) be a cut-off 
function such that = 1 in Q{zo,R) and |Vv9| < C/i?, \dtip\ < C/R^ in 
(5(-2o, 2-R). Choose v = v^^U;^ as a test function in ( l35i) . then 

+2 / fi{9^)(pDu^ : {u^ (g) Vv?) =1+1 v^^f • u,, + 2 / Pi.v'Uj, ■ Vy? 



where / corresponds to the convective term 

1=1 A^^(u,,) (g) : (v9^Vu^ + 2v?Vv? ® Uj,)(i2; 
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= / (MuM ■ Vy?) \u^\'^^dz. 

JQt 

Applying the assumption ([5]) and Holder and Young inequalities, it follows 

1 ]-{ip^\\i^\^){t)dx + I if^\D\ii,f'dxdT < I ip\\i^\^\dtip\ + 
JQ 2 Jq^ Jq^ 

+6 [ (^|Vu,p + C(5,/i#) / ip\u,\^\V^\^ + \I\ + 

jQt JQt 

+ \l ^'If|' + ^/ + |p,|("+2)/n^(«+2)/(2n) + 

2 JQt ^JQt ^ + 2 jQt 

I ^ ^ /" ,o("+2)/4|„ |(n+2)/2|^ |(n+2)/2 

By Korn inequality the following estimate holds 

(^2|Vu,|2<2/ ip^\Du,\^ + A[ |V<^|V.p. 



According to the properties of (f it arises 



JQ{zo,R) JQ{za,2R) ^ JQ{zq,2R) 

+ |/| + C / |f prfz + ^ / 

JQ(zn,2R) n + 2 jQ(zn,2K) 



lQizo,2R) If-^ ^ JQ{z„,2R) 

observing that for _R < 1 we have R"^ > = i?'^/2+(n+2)/2_ gj^ce 



|/|<^ / |A^,(u,)||u,|2dz< 
-n^ JQ(zn,2R) 



lQ(zo,2R) 

c c 

then we conclude (1381). 



3.5 Existence of weak solutions 

From estimates dH]), (EDD-lEI]), (I33D-(I31D and (ETj), independent on z/, we can 
extract a subsequence, still denoted by (yiy.py.Oy), verifying fl3T]) - fl32l) and 

^ u weakly* in L°°(0, T; V{Sl))] 
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Vu 


weakly in 




(39) 


dtn^ 


dtu 


weakly in 






- 


-> u 


strongly in 


V\Qt), for 1 < m < 2(n + 2)/n; 


(40) 




^ e 


weakly in 


L^(0, T- Wl'\n)), for 1< g < 2 - 


n/{n + 1); 


eu 




strongly in 


L'^iQr), for 1 < m < q{n + l)/n; 


(41) 


Pv 




weakly in 






Pu,l - 


^pi 


strongly in 


L™(gT), for 1 < m < (n + 2)/n; 


(42) 


Pu,2 - 


^P2 


weakly in 


L\Qt). 





Since div u^, = and u^\dnx]o,T[ = 0, it follows from the definition and 
properties of Helmoltz decomposition that 

G L°°(0, T; L"'{n)) for all m G [1, oo), 
Muiu^) u strongly in L"'(Qt), for 1 < m < 2{n + 2)/n. 

Then, the initial value problem (!3T]) - (!32|) passes to the limit as u tends 
to infinity, concluding the required problem ©-(jH]). 

3.6 Local energy inequalities 

3.6.1 Proof of the local energy inequality ([9]) 

Let ip G C^{Qt) and choose v = (p'^{ui, — a) as a test function in (|35|) for an 
arbitrary a G M", arguing as in Section [231 then 



2 

|2 



^y{u^-a.)\\l^{t)+ I fi{9^)<^'^\Du^\'^dxdT 



(p\u^ — a| dt^pdxdr — / fi{6^)ipDu^ : (u,^ — a) (g) Vifdxdr + 

JQt 

+ / (Muiu^) ■ Vv?) u,, ■ (u,, - SL)ipdxdr 
jQt 

+ / V9^f ■ (u^ - a) + 2 / Pu(~p{viu - a) ■ Vipdxdr. 

jQt JQt 

Next, arguing as in Section 13.51 we can pass to the limit as u tends to 
infinity concluding 
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3.6.2 Proof of the local energy inequality ( llOj) 

Going back to the solution (u^, p^, 6^) G (L~(0, T; L\n))r]L^{0, T; Wl'\n)))x 
L^{0,T;W^'^{n)) X £ of the initial value problem ([IH])-([l9D we can clioose 
= e^iC + (^^)2)-i/2^ e L2(0,T; Wo'\n))), for C > O and a non-negative 
function ip G C^{Qt), as a test function in (HM . First, we calculate sepa- 
rately the following terms 



(9*^^,0)= / (y/(+{e^y^) {x,t)dx- / ^/cTW?dtijix,T)dxdT 



^/CTW^M,{u^) ■ Vijdxdr 

jQt 

(V^^, V0) = C J^^^QNYyri^dxdT + V [VC + W?) ■ VijdxdT 
= C J^^j0N\^2^dxdT - y^CTWfA^dxdr. 

Thus, we conclude 



= / + {d''? {dt^ + kAij + Muiu^) ■ V^) dxdr + 

JQt 

+ f fi(9^)\Du''\^^l^=dxdT. (43) 

Since ip ^ C^{Qt) such that -0 > we can pass successively to the limit 
as N tends to infinity and e and u tend to zero, obtaining fllOj) . Indeed, 
still denoting by (u^,p^,6'^) the subsequence extracted from the solutions 
(u^,p^,^^) of (HSD-dED it verifies 

Ue weakly* in L~(0, T; L'^{n)); 

Vu^ ^ Vug weakly in Li'^(Qt)'-i 

dtu^ dtu, weakly in L'^{0, T; W'^'^in)); 

u, strongly in L'^{Qt), for 1 < m < 2(n + 2)/n; 
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e, weakly in ^"(0, T; W^o'^(n)), for 1< g < 2 - n/{n + 1); 
6^ 9^ strongly in L'^{Qt), for 1 < m < q{n + l)/n; 
Pe weakly in L^(0, T; ^^^•^(^^)); 



where {u^,Pe, O^,) is a solution of (BBl) - (|7n) . Moreover, considering that ^ > 0, 

(CT(n¥^ ^ (CT#)^ w^-kly in L^(Q^), 
the lower semicontinuity property of the L^-norm and 

fjN a 

/i(g^) , = ^ ^Ji{e,)^^= *-weakly in "L^^iQr), 

Wc + Wf Wc + W 

and (!25|) hold, then the equality (143|) delivers to the local energy inequality 
([To]) satisfied by (vie.Pe.Oe). 

Next, still denoting by (u£,p£,6'e) the subsequence extracted from the 
solutions (u£,Pe,6'e) of (|2B]) - (I27|) it verifies 

u, weakly* in L~(0, T; L2(fi)); 

Ue My weakly in ^^(0, T; Wo'^(^])); 

dtVL, dtVLy weakly in ^^(0, T; W-^'2(fi)); 

Ug ^ \iy strongly in L'"(Qt), for 1 < m < 2(r?, + 2)/?7,; 

Oe By weakly in L'?(0, T; IVo'^(n)), for 1< g < 2 - n/(n + 1); 

9e 9u strongly in L"^{Qt), for 1 < m < q{n + l)/n; 

^ Pi, weakly in L^{Qt)- 

where {ui,,p^,9i,) is a solution of fl3Tl) - fl32l) . Analogously to the above argu- 
ment the local energy inequality ffTOl) arises for (ui,,pi,,9i,). 

In conclusion, the convergences in Section 13.51 imply that the limit 9 
verifies ^ > in Qt and the local energy inequality flTUl) holds for the weak 
solution {u,p,9) such that verifies 

3.6.3 Proof of the local energy inequality ( llip 

Proceeding as in Section 13.6.21 we go back to the existence of the solution 
(u^,p^,^^) G {L^{0,T;L^{n))nL\0,T;Wl'\n))) X L^{0,T;W^'\n)) X S 
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of the initial value problem fll8p - fll9p . Next we can choose (p = (1+6^) & 
L^iO,T; Wo'^{n))), for < ^ < 1 and a non-negative function ^ G C^iQr), 
as a test function in f|T9l) . Now, calculating separately the following terms 



1 - 






1 


1 






1 


1 




-e 


/ 

JQt 



(9,0^,0) = I {l+9'')'-^^{x,t)dx- 

'l+e^y-^^t^{x,T)dxdT 



Thus, we conclude 



--^ / il + e^Y-^{dtij + kAi} + Muiu^)-V^)dxdT. 

^~ ^ JQt 

Since ■?/' G C^{Qt) such that > we can pass successively to the limit as 
tends to infinity and e and u tend to zero, obtaining (fTTj) . 



4 Regularity of 9 (Proposition 12.11) 



Let (u, p, 9) be a weak solution in accordance to Theorem I2.H that is, it 
satisfies in the sense of distributions 

dtO - kAe = -u-Ve + |Dup in Qt. (44) 



Thanks to Theorem 12.11 we have 



iVu^ G L^+^(gT). 
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Thus we get 



1 n 1 

X 4(n + 2)(l + ej q 



Since g < (n + 2)/(n + 1) it follows 

4(n + 2)(l + e) 

1 < X < ^ —, for every e > 0. 

5n + 4 + 4e(n + l) 

In particular, x>l + eife<(4 — n)/ (4n + 4). 
Let us split the proof into two cases. 

n—2 Since we have u-VO G L^^^'^^Qt) supposing e < 1/6, the classical theory 
for the heat equation leads 9 G W^i+^e,ioc(<5T)- 

n=3 The classical theory for the heat equation (jUj) leads 6 G W^^'ioc(Qt) *■ 
^5../{5-.<)(0^r; Vr^'5V(5— )(^)). Thus we get u ■ V0 G /^io+'I^t) sup- 
posing e < 1/6, and the bootstrap argument guarantees the required 
result. 

Finally, the last assertion is valid from the embedding 

References 

[1] A. Arkhipova and O. Ladyzhenskaya, On inhomogeneous incompress- 
ible fluids and reverse Holder inequalities. Annali della Scuola Normale 
Superiore di Pisa - Classe di Scienze, Ser. 4, 25: 1-2 (1997), 51-67. 

[2] A. Arkhipova and O. Ladyzhenskaya, On a modification of Gehring's 
Lemma. Zapiski Nauchnykh Seminarov POMI 259 (1999), 7-18; English 
transl.: J. Math. Sci. (N.Y.) 109 :5 (2002), 1805-1813. 

[3] L. Boccardo, A. Dall'aglio, T. Gallouet and L. Orsina, Nonlinear 
parabolic equations with measure data. Journal of Functional Analy- 
sis 147 (1997), 237-258. 

[4] L. Boccardo and T. Gallouet, Non-linear elliptic and parabolic equations 
involving measure. Journal of Functional Analysis 87 (1989), 149-169. 



18 



[5] G. Boling and Y. Guangwei, On the suitable weak solutions for the 
Cauchy problem of the Boussinesq equations. Nonlinear Analysis: The- 
ory, Methods & Applications 26 :8 (1996), 1367-1385. 

[6] M. Buhcek, L. Consiglieri and J. Malek, Slip boundary effects on 
unsteady flows of incompressible viscous heat conducting fluids with 
a non-linear internal energy-temperature relationship. Preprint CMA- 
FUL, Pre-2007-014. 

[7] L. Caffarelli, R. Kohn and L. Niremberg, Partial regularity of suitable 
weak solutions of the Navier-Stokes equations. Comm. Pure Appl. Math. 
35 :6 (1982), 771-831. 

[8] L. Gonsiglieri, Weak solutions for a class of non-Newtonian fluids with 
energy transfer. J. Math. Fluid Mechanics 2 (2000), 267-293. 

[9] L. Consiglieri, Friction boundary conditions on thermal incompressible 
viscous flows. Annali di Matematica Pura ed Applicata 187 :4 (2008), 
647-665. 

[10] L. Consiglieri, Regularity for the Navier-Stokes-Fourier system. Differ- 
ential Equations & Applications 1 :4 (2009), 583-604. 

[11] L. Consiglieri and T. Shilkin, Regularity to stationary weak solutions 
for generalized Newtonian fluids with energy transfer. Zapiski Nauchnyh 
Semmarov POMI 271 (2000), 122-150; English transl.: J. Math. Sci. 
(N.Y.) 115 (2003), 2771-2788. 

[12] D. Donatelli, On the artificial compressibility method for the Navier 
Stokes Fourier system. arXiv:0807.3842vl [math.AP] 

[13] G.P. Galdi, An introduction to the mathematical theory of the Navier- 
Stokes equations. Linearized steady problems. Springer Tracts in Natural 
Philosophy 38, New York 1994. 

[14] F.W. Gehring, The L^-integrability of the partial derivatives of a quasi 
conformal mappings. Acta Math. 130 (1973), 265-277. 

[15] J.-L. Guermond, Faedo-Galerkin weak solutions of the Navier-Stokes 
equations with Dirichlct boundary conditions are suitable. J. Math. 
Pures Appl. 88 (2007), 87-106. 



19 



[16] S. Gustafson, K. Kang and T.-P. Tsai, Regularity criteria for suitable 
weak solutions of the Navier- Stokes equations near the boundary. J. 
Differential Equations 226 (2006), 594-618. 

[17] T. Hishida, Existence and regularizing properties of solutions for the 
nonstationary convection problem. Funkcialaj Ekvacioj 34 (1991), 449- 
474. 

[18] O.A. Ladyzhenskaya, Mathematical problems in the dynamics of a vis- 
cous incompressible fluid. 2nd rev. aug. cd., "Nauka", Moscow, 1970: 
English transl. of 1st ed.,The mathematical theory of viscous incom- 
pressible flow. Gordon and Breach, New York 1969. 

[19] O.A. Ladyzhenskaya, V.A. Solonnikov and N.N. Uraltseva, Linear and 
quasilinear equations of parabolic type. Translations of Mathematical 
Monographs, Vol. 23 American Mathematical Society, Providence, R.L 
1967. 

[20] F. Lin, A new proof of the Caffarelli-Kohn-Niremberg theorem. Comm. 
Pure Appl. Math. 51 :3 (1998), 241-257. 

[21] J.L. Lions, Quelques methodes de resolution des problemes aux limites 
non lineaires. Dunod ct Gauthier-Villars, Paris 1969. 

[22] J. Malek, J. Necas, M. Rokyta and M. Ruzicka, Weak and Measure- 
valued solutions to evolutionary PDEs. Chapman and Hall, London 
1996. 

[23] J. Naumann, On the existence of weak solutions to the equations of 
non-stationary motion of heat-conducting incompressible viscous fluids. 
Math. Meth. Appl. Sci. 29 (2006), 1883-1906. 

[24] J. Naumann and J. Wolf, On the interior regularity of weak solutions to 
the non-stationary Stokes system. J. Glob. Optim. 40 (2008), 277-288. 

[25] T. Roubi'cek, On non-Newtonian fluids with energy transfer. 
J. Math. Fluid Mech. 11 (2009), 110-125 (published on line in vol.9, 
no.l (June 2007), 1-16.). 

[26] V. Scheffer, Hausdorff measure and the Navier-Stokes equations. Comm. 
Math. Phys. 55 :2 (1977), 97-112. 



20 



[27] V. Scheffer, The Navicr-Stokes equations on a bounded domain. Comm. 
Math. Phys. 73 (1980), 1-42. 

[28] E.W. Stredulinsky, Higher integrabihty from reverse Holder inequahties. 
Indiana Univ. Math. J. 29 (1980), 408-417. 

[29] J. Simon, Compact sets in the space. Annali Mat. Pura Appl.IV 146 
(1987), 65-96. 

[30] W. Walter, Ordinary differential equations. Springer- Verlag, New York, 
Inc. 1998. 



21 



